Flow Induced Organization and Memory of a Vortex Lattice 
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We report on experiments probing the evolution of a vortex state in response to a driving current in 
2H-NbSe2 crystals. By following the vortex motion with fast transport measurements we find that 
the current enables the system to reorganize and access new configurations. During this process 
the system exhibits a long-term memory: if the current is turned off the vortices freeze in place 
remembering their prior motion. When the current is restored the motion resumes where it stopped. 
The experiments provide evidence for a dynamically driven structural change of the vortex lattice 
and a corresponding dynamic phase diagram that contains a previously unknown regime where the 
critical current can be either increased or decreased by applying an appropriate driving current. 



PACS numbers: 74.60. Ge, 74.60. Jg, 74.60Ec 

Vortices in a type II superconductor exhibit remark- 
ably rich and sometimes surprising coUective behavior. 
The underlying physics is a delicate balance between two 
competing forces: one is the mutual interaction between 
vortices which promotes order and the other is pinning, 
caused by randomly distributed material defects, which 
favors disorder . A striking example of this competi- 
tion is the " peak effect" - a sudden increase in the critical 
current just below the superconducting transition tem- 
perature, Tc, attributed to a violent structural reorga- 
nization when pinning overcomes the ordering influence 
of interactions . Another outcome of this balance is 
that the vortices are subject to an energy landscape with 
many metastable states leading to history dependence 
and nonlinear behavior [p[-p^. For example a vortex ar- 
ray prepared by field cooling, where the magnetic field is 
applied before cooling through Tc , becomes trapped in a 
metastable state that has a higher critical current and is 
more disordered than one prepared by zero field cooling 
- where the field is applied after cooling |^,^. Setting 
the vortices in motion with an external current releases 
them to explore the energy landscape and to access new 
states , but thus far the dynamics of this process was 
inaccessible to experiments. 

To probe the vortex dynamics we have devised a tech- 
nique that uses fast current ramps and pulses to track 
the evolution of a vortex lattice as it unpins and starts 
moving. We find that when the vortices are set in motion 
they organize into a new state whose properties are de- 
termined by the driving current. The experiments lead 
to a dynamic phase diagram containing four regimes as 
a function of increasing current: A) pinned vortex array 
at the lowest currents; B) at higher currents a region of 
motional reorganization and repinning where the vortices 
settle into a state whose critical current is equal to the 
driving current; C) at still higher currents channels of or- 
dered moving vortices cut through domains of disordered 
stationary ones; D) at the highest currents all the vor- 
tices unpin and form an ordered state with low critical 



current. 

The vortex motion is detected by measuring the lon- 
gitudinal voltage across the sample V— vBl, where v is 
the vortex velocity, B the magnetic field and I the dis- 
tance between the voltage contacts. As long as the vor- 
tex array is pinned the current flow is dissipationless and 
only when a critical current Ic is exceeded, unpinning it, 
will a voltage appear. It was first shown by Larkin and 
Ovchinnikov |^ that for weak pinning superconductors Ic 
is a measure of the degree of order in the vortex system. 
In this case, as was confirmed by recent SANS and 
transport measurements [||, the critical current density 
is inversely proportional to the size of an ordered vortex 
domain and Ic can be used, as was done here, to quantify 
the degree of vortex order. 

The data were acquired on a Fe (200ppm) doped single 
crystal of 2H-NbSe2 of size of (3x1x0.015) mm^ and Tc 
— 5.95K. Our measurements employed a standard four 
probe technique with low resistance contacts made with 
Ago. lino. 9 solder. The current ramps and pulses were 
generated with a Wavetek 164 and the voltage response 
was measured with a fast (INA103) amplifier. The mag- 
netic field was kept along the c axis of the sample and the 
current in the a-b plane. All the data were taken after 
the sample is zero field cooled (ZFC). Previous experi- 
ments on similar samples revealed novel types of nonlin- 
ear behavior including frequency memory and long glassy 
relaxation times that were observed in the lower part of 
the peak effect region p^ . In the present experiments 
we find that this is also the region where the value of the 
critical current changes during the measurement. 

Fig.l shows the response to two types of current ramps, 
fast (FCR -ramping rate of 200A/s), and slow current 
ramps (SCR -ramping rate of ImA/s). The critical cur- 
rent was defined as the current at which the voltage 
reaches 1/iV. The inset in Fig.l, which shows the current- 
voltage {I ~V) characteristics for both ramps, illustrates 
that the response depends on the speed of measurement. 
In particular the critical current Icf, obtained with the 



1 



120 




4.2 4.4 4.6 4.8 5.0 5.2 



T(K) 

FIG. 1. Temperature dependence of Ic detected by slow 
(SCR) and fast (FCR) current ramps. Inset: I — V curves 
in response to SCR (ImA/s) and FCR (200A/s). The labels 
A,B,C,D illustrate the four regions of the dynamic phase di- 
agram. 

FCR, is lower than Ics , the value obtained for the SCR. 
The temperature dependence of Ics and Icf is shown in 
the main panel of Fig.l. A strong peak effect as re- 
ported in Ref.[5] is seen in both data sets. The fact that 
Icf < Ics, for all data points means that the ZFC state 
probed by the FCR is always more ordered than the SCR 
state, leading to the conclusion that during the SCR the 
vortex array has time to reorganize into a more strongly 
pinned state. This organization does not however pro- 
duce a voltage signal in the SCR measurements because 
it decays too fast to be resolved. 

To observe the reorganization directly we carried out 
fast measurements with 2/isec temporal resolution, which 
enabled us to follow the decay of the voltage response to 
a current step / > Icf as shown in Fig. 2. We note that 
upon applying the current step the vortex response is in- 
stantaneous as expected of essentially massless objects 
[ p^ when they respond elastically. As time progresses 
the voltage decreases despite the fact that the driving 
current is kept constant. After a sufficiently long waiting 
time the voltage decays to zero for driving currents in 
the range Icf < I < Ics- The decay is also observed for 
driving currents in excess of Ics indicating that the mov- 
ing vortex array organizes itself into less ordered state 
even though it is not coming to a full stop. In either case 
the organization proceeds in two steps. On long time 
scales there is a clear decay which exhibits a stretched 
exponential time dependence exp(-t/r)" where both the 
characteristic time scale, r and the exponent a depend on 
current and temperature. This type of relaxation is typ- 
ical of the dynamics observed in glasses. By contrast, on 
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FIG. 2. Time dependence of response to a current step. 
The inset shows the response to steps of different amplitudes. 

short time scales the response is almost constant showing 
no measurable decay. This is consistent with the obser- 
vation of Ref. [8] that the vortex pattern does not change 
instantaneously. It is this delay before the decay sets in 
that makes it possible to observe the pristine vortex state 
- before it has time to reorder - by using a fast current 
ramp technique. 

During the reorganization a striking long-term mem- 
ory is observed. A typical result is shown in Fig. 3 (a) 
where we plot the response to a current pulse of am- 
plitude /. When the current is interrupted during the 
decay, the voltage drops to zero instantaneously. Sub- 
sequently, when the current is restored after a waiting 
time with zero current, the voltage decay resumes where 
it stopped the moment the current was turned off. The 
waiting times for which the lattice remembers its state 
of motion prior to being stopped can be very long (we 
checked waiting times of up to 20 hours) provided the 
field and temperature remain unchanged. Other effects 
associated with quenching of the vortex motion where 
observed in experiments on the structure 0JlJ,|l^ and 
transport properties pof . 

Long-term memory occurs in systems that can be 
trapped indefinitely in a metastable state. In the case 
of a vortex array this implies that the energy scale of 
the pinning potential is much larger than that of thermal 
fluctuations. The current lowers the barriers of the pin- 
ning potential and sets the vortices in motion allowing 
them to explore new states and eventually settle into a 
more strongly pinned state. It is thus the applied cur- 
rent rather than the temperature that drives the system 
towards equilibrium. How robust are these metastable 
states or equivalently how large can the current be before 
the vortices become untrapped? To explore the interplay 
between current and pinning we carried out the experi- 
ment shown in Fig. 3(b). The ZFC vortex lattice is driven 
with a current pulse of amplitude / for 1.6nisec during 
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FIG. 3. Memory of a vortex array, (a) During the decay 
the current is removed and later it is restored to its initial 
value. The lines and the open circles are the current and 
the response respectively, (b) SV as a function of Iq. The 
definitions of SV and Iq are shown in inset. 



which the response decays correspondingly. The current 
is then suddenly reduced to a lower value Iq for 6msec 
during which the voltage decays at a different rate, and 
subsequently it is restored to its original value. Plotting 
SY, the voltage decay during the reduced current pulse, 
as a function of Jq , we find that there is indeed a finite 
value of Iq below which the system does not evolve with 
time and retains its memory. This value of the current 
coincides with the critical current of the vortex lattice at 
the moment at which the reduced current is applied. 

The memory can be a powerful tool in the study of 
vortex dynamics. Most probes of the vortex structure re- 
quire long times to produce one data point or a picture, 
typically 45 seconds in scanning Hall probe microscopy 
PH and 15 minutes in SANS 0. It would be impossible 
to observe the motional reorganization with these tech- 
niques. But as our experiments have shown, the vortex 
structure can be frozen at any moment during its motion 
by removing the driving current. The quenched vortex 
system maintains its structure and can then be probed 
directly. We used this method here to observe the vortex 




10 15 20 25 30 35 
Current (mA) 



40 
30 

< 
E 

~„ 20 
10 



T= 4.80K 



B= 0.4T 



<>oOo- 

% 20mA 



,4* □□ ^ □ 



A- B A' — C 



40 
30 

<■ 
E 

20:1^ 
10 



1 0"''l O'^l 0"^ 0"M 0° 1 0' 10 20 30 40 50 
^0 Driving current (mA) 

FIG. 4. (a) I — V curves obtained with a fast current 
ramp following decay periods to=0, 3ms and 3s with current 
J=33mA. The definition of to is shown in the inset, (b) Evo- 
lution of critical current with time in the presence of driving 
currents / = 20, 33, and 52mA. (c) Data for 7c obtained at 
to = 3s and various driving currents. The symbols A, B, C 
represent the three regimes of the dynamic phase diagram. 
The dashed line is a guide to the eye. 



reorganization by quenching the system at various times 
during its motion and measuring its critical current with 
the PGR method. 

As an example, we show in Pig. 4(a) the response to 
an PGR at decay times tg — 0, 3ms and 3s. Glearly, 
Ic increases with increasing decay time. In Pig. 4(b) we 
show the dependence of Ic on decay time for driving 
currents / = 20, 33 and 52mA. Por current amplitudes 
Icf < I ^ Ics, the critical current increases with decay 
time and saturates at a value which is approximately 
equal to the driving current. In Pig. 4(c) we show the 
dependence of Ic on / following a long decay time of tg 
= 3s. As long as Icf < I < Ics we find that Ic ~ /, in 
other words the system can be prepared with any de- 
sired critical current in this range simply by priming it 
with the appropriate driving current. Por I > Ics, there 
is an abrupt drop in Ic to a value only slightly above that 
in the ordered state. We interpret this drop in terms of 
the formation of a channel of unpinned and motionallay 
ordered vortices. The fact that critical current is slightly 
higher than that of the ZPG state may indicate that, as 
predicted by theory [|lO|,0 , the motionally ordered state 
resembles a smectic rather than a crystal. 
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We now show that for / > Ics our data is consistent 
with channel formation and estimate the number of vor- 
tices in channels. First we interpret the information con- 
tained in the I — V curves in terms of vortex motion: 
Ic gives the force required to set the vortices in motion 
and the cordal resistance, defined as Rc ~ V/{I — Ic), 
measures the number of moving vortices in response to 
a current /. When all the vortices are moving Rc takes 
the maximum value, given by the Bardeen-Stephen resis- 
tance Rbs ~ Rni'^o/ Hc2)'^, where i?„ is the normal 
state resistance, $o the fundamental unit of flux, Hc2 the 
upper critical field and n is the vortex density. If only 
a fraction of the vortices move, Rc is reduced propor- 
tionately. We note that for the FCR data in the inset of 
Fig.l the cordal resistance is close to the maximum value, 
Rl= 8.75ml7 = QMRbs which means that essentially all 
the vortices move. By contrast the I — V curve obtained 
with the SCR for the same initial vortex state is different: 
it has a much higher critical current and once the vor- 
tex motion is detectable the slope increases with current. 
In this case the number of moving vortices depends on 
the current and can be estimated by using the result of 
Fig. 4(c) that for I > Ics , the moving vortex array is in a 
state with critical current ^ Icf- Therefore, the ratio of 
voltage response obtained with the SCR to that obtained 
with the FCR Vs{I)/Vf{I) = Rl/R{ is equal to the ra- 
tio of moving vortices in the two measurements. Since in 
the FCR case the entire array is moving Vs{I) / Vf{I) is 
in fact a direct measure the fraction of moving vortices. 
This result leads to the conclusion that for Ics < I < ZIcs 
not all the vortices are participating in the motion and 
those that do move form channels in which the vortices 
are ordered (because the critical current is close to Icf 
), consistent with our assumption. The fraction of mov- 
ing vortices in the SCR increases gradually from zero at 
/ = Ics to 1 when the channels have engulfed the entire 
sample forming an ordered state. 

Based on these results we identify four regions in the 
dynamic phase diagram as illustrated in Fig.l. In region 
A where I < Icf , the vortex lattice stays pinned and 
ordered. Region B, where Icf < I < Ics, is the regime 
of motional reorganization and repinning. Here the re- 
sponse decays with time as parts of the lattice break loose 
settling into stronger pinned configurations until the mo- 
tion comes to a halt when the vortices reach a state with 
Ic equal to the driving current. This state contains at 
least one contiguous region of disordered vortices that 
provides a low resistance current path along the sam- 
ple. Region C, where Ics < I < 2Ics , is the plastic flow 
regime in which channels of moving ordered vortices cut- 
ting across the sample are embedded in domains of dis- 
ordered stationary vortices. In region D, / > 2Ics , the 
entire vortex array unpins and becomes motionally or- 
dered. 

We conclude that the current assists the system in 
choosing a configuration in the energy landscape formed 



by the competition between interactions and pinning. 
The current assisted rearrangement can also account for 
the fact that the ZFC vortex state has the lowest critical 
current: as the field is ramped up in the superconduct- 
ing state, it is initially expelled from the interior of the 
sample by the screening currents. The vortices can pen- 
etrate only when the screening current density exceeds 
the critical value, and therefore the configuration with 
the lowest critical current- the ordered state -will 
be the first to penetrate the sample. For the experi- 
ments discussed here the typical ordered domain in the 
ZFC state contains ~ 10'^ vortices. Once the vortices are 
inside the sample an externally applied current assists 
them in finding a lower energy configuration by adapting 
to the pinning potential and becoming more disordered. 
With increasing current, as the pinning barriers are ef- 
fectively lowered, larger portions of the pinscape can be 
explored resulting in progressively more strongly pinned 
states. Until at / > Ics a channel forms across the sam- 
ple in which no barriers are encountered allowing free 
passage of moving, ordered vortices. The motional reor- 
ganization does not occur immediately after the array is 
set in motion but after a delay during which it maintains 
its prior configuration. This initial hesitation together 
with the long term memory makes it possible to directly 
probe these processes with transport measurements and 
in principle with other techniques as well. 
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